The Travelling Cluster Approximation for Strong Correlation Models of Lattice 

Fermions Coupled to Classical Fields 
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We suggest and implement a new Monte Carlo strategy for correlated models involving fermions 
strongly coupled to classical degrees of freedom, with accurate handling of quenched disorder as 
well. Current methods iteratively diagonalise the full Hamiltonian for a system of A'^ sites with 
computation time rjv ~ A^'*. This limits achievable sizes to TV ~ 100. In our method the energy 
cost of a Monte Carlo update is computed from the Hamiltonian of a cluster, of size Nc, constructed 
around the reference site, and embedded in the larger system. As MC steps sweep over the system, 
the cluster Hamiltonian also moves, being reconstructed at each site where an update is attempted. 
In this method TN,Na ~ NN^. Our results are obviously exact when A'^c = A'^, and converge quickly 
to this asymptote with increasing Nc . The accuracy improves in systems where the effective disorder 
seen by the fermions is large. We provide results of preliminary calculations on the Holstein model 
and the Double Exchange model. The 'locality' of the energy cost, as evidenced by our results, 
suggests that several important but inaccessible problems can now be handled with control. 



I. INTRODUCTION 



The equilibrium physics of classical interacting systems 
is by now very well understood. Quantum many body 
problems, involving strong interactions, however, remain 
difBcult to solve with control. The focus of strong cor- 
relation theory is on devising methods to handle these 
problems. Techniques like density matrix renormalisa- 
tion group^ (DMRG) and dynamical mean field theory^ 
(DMFT), for example, represent an advance in this direc- 
tion. Problems where fermions are coupled to 'classical 
fields', e.g, large S spins, or 'adiabatic' phonons, are at 
an intermediate level of difficulty between purely clas- 
sical systems and quantum many body problems. The 
quantum degrees of freedom are not directly interacting, 
so the difficulty with an exponentially growing Hilbert 
space is absent, but "annealing" the classical variables is 
much more difficult compared to purely classical systems. 

The adiabatic approximation, whereby some degrees 
of freedom are treated as classical, is not novel. Several 
problems have been solved in the past by making this 
approximation, e.g, in electron-phonon systems^, or, in 
a different context, in the Car-Parrinello method^, han- 
dling coupled electronic and ionic degrees of freedom. 
The recent interest lies in the application of this approach 
to several strong coupling lattice fermion models, and 
some degree of success in understanding complex ma- 
terials. Millis and coworkers^ studied electrons coupled 
to classical spin and lattice ('phonon') degrees of free- 
dom, to construct an initial theory of the manganites 
using DMFT. The approach was taken much farther by 
Dagotto and coworkers^ using 'real space' Monte Carlo 
techniques to study ordering phenomena, phase coexis- 
tence, and disorder effects in a large family of correlation 
models pertinent to the manganites. The method has 



been used extensively also to explore magnetism in dou- 
ble exchange (DE) based models^. For diluted magnetic 
semiconductors (DMS) too much of the physics has been 
clarified by methods which treat the doped magnetic 
moment as classical*. We ourselves have used the ap- 
proach to study magnetism, insulator-metal transitions 
and nanoscale phase coexistence in disordered correlated 
electron models^'^". 

The adiabatic limit simplifies the many body problem 
by casting it in the form of "non interacting" fermions 
in the background of classical variables, {x} say, but de- 
termining the distribution P{x} involves an expensive 
computation. In the absence of any small parameter to 
simplify the problem, computing P{x} requires iterative 
diagonalisation^^ of the fermion Hamiltonian and, for an 
N site system, the computation time, tjv, increases as 
N'^. We will describe the standard exact diagonalisa- 
tion based Monte Carlo (ED-MC) in the next section, 
here we just note that the accessible sizes, N ~ 100, 
within ED-MC, severely limits the ability of the method 
to resolve the outstanding issues relating to transport, 
metal-insulator transitions, and the effect of disorder in 
correlated systems. 

There have been some attempts at overcoming the se- 
vere finite size constraint in ED-MC. (i) Instead of exact 
diagonalisation, it has been proposed that the energy of 
fermions in the classical background can be estimated 
by moment expansion of the density of states. This, in 
principle, allows access to iV ^ 10'^, and has been used 
to study the clean^^ and disordered^'^ DE model, {ii) A 
'hybrid' Monte Carlo method, using dynamical evolution 
of the classical variables, has been tried out^^ for a model 
of competing DE and superexchange. (iii) We have pro- 
posed a scheme^^, in the context of double exchange, 
where the energy associated with the spin configuration 
can be approximated by an explicit classical Hamiltonian 



with couplings determined from a solution of the fermion 
problem. 

While the approximations above have allowed some 
advance in the context of double exchange, there is no 
equivalent method available for handling phonon degrees 
of freedom, or the combination of phonons and spins 
(as relevant to manganites), or dilute strong coupling 
systems like the DMS. There is the need for a general 
and computationally transparent method that can handle 
models with arbitrary coupling and disorder, and sys- 
tematically approach the 'exact' answer. This paper pro- 
poses such a scheme. We employ a variant of the exact 
diagonalisation strategy using an embedded (travelling) 
cluster, that estimates the energy cost of a Monte Carlo 
move by diagonalising the smaller cluster rather than the 
full Hamiltonian. Since size limitations are most severe 
in three dimension {3d), and it is physically also the most 
relevant, we benchmark our method directly in 3d, where 
the test is most stringent. 

We will study two models, in 3d, to provide some per- 
formance benchmarks on the travelling cluster approxi- 
mation (TCA). These are (i) the (disordered) Holstcin 
model, and (m) the double exchange model. They are, 
respectively: 

i^i = -t ^ cjcj + X!*^*'* ^ ^^'^^ ^ ^ X! "^^^^ '^ ^^ 
H2^-tY^ cl^Cj^ + ^(e, - n)n^ ~ Jh^ S,.o-j (1) 
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The t are nearest neighbour hopping on a three di- 
mensional lattice. In Hi, the on site binary disorder 
Ci assumes value ±A, /x is the chemical potential and 
Ui = cjci is the electron density operator (for spin- 
less fcrmions), coupling to the phonon coordinate Xi. 
Hk =■ {K/2) ^^ xf where if = 1 is the stiffness of the 
phononic oscillators. In H2, Jh is the Hunds coupling. 
The model is defined with 'spinfuU' fcrmions, but since 
we will use Jh /t ^ 00 it will also lead to an effective 
spinless fermion problem. We set t = 1, fixing our basic 
energy scale, and also h— I. 



II. METHOD 

Let us start with the T = case to clarify the usual ap- 
proach to these problems. Since we have earlier discussed 
the ED-MC method in detail^^ we only provide a brief 
outline here. There are several applications of ED-MC 
in the context of manganite related models^^. 

There are two (related) difficulties in solving strong 
coupling adiabatic problems. {i) The probability of 
occurence of a classical configuration is not explicitly 
known, and is governed by the fermion free energy. Gen- 
erating these configurations involves the N^ cost spec- 
ified earlier. Let us call this the "annealing problem". 
{ii) Even for a specified classical configuration, obtained 



via some annealing technique, the electronic properties, 
e.g, the resistivity, involves computing fermionic correla- 
tion functions in a non trivial background. Since there 
is no analytic theory for non interacting fcrmions in an 
arbitrary 'landscape', transport calculations have to im- 
plement the Kubo formulae exactly. Our innovation in 
this paper is on the annealing problem, we still depend on 
a numerical implementation of linear response theory^^ 
to solve the transport problem on large lattices. 

Let us start with the annealing problem. If we are 
at T = 0, the background in which the fermions move 
can be determined by minimising the total energy {H) 
with respect to the classical variables. Denote the classi- 
cal configuration as {??i,?72, ..} = {r/}, where rji ~ ?7(Ri) 
and, in case of multiple classical variables at each site, rji 
represents the full set of variables Xi, Si, etc, at that site. 
The key task is to determine {H) — £{rj}. 

If the coupling between the classical and quantum vari- 
ables is large, there is no perturbative result for the 
fermion energy £{t]}, and therefore no explicit functional 
that we can minimise. This is where MC is used. The 
exact diagonalisation based method uses the following 
strategy: («) Set up an arbitrary configuration {77} and 
compute the energy £{vi\. The fermion contribution is 
estimated by direct diagonalisation while the classical 
contribution, Kxf say, is explicit, {ii) Attempt an up- 
date, say at site R^, by changing rji -^ ri[. Compute the 
energy £{ri'}. {Hi) If Af ~ ^iv'} ^ ^iv} < 0: accept 
the move, if Af > 0, and T 7^ 0, accept the move with 
probability oc e^ '"^. {iv) Sweep over the system, ini- 
tially to reach equilibrium and then to compute thermal 
averages. 

The method above is simply a use of the Metropolis 
algorithm, with the complication of an expensive diago- 
nalisation for every microscopic update. Since each local 
update involves computational effort ^ N^ , the cost of 
sweeping over the system leads to tjv ~ N^. One has 
to multiply this with the cost of thermal averaging, and 
disorder average (if needed) . 

We were motivated to ask if it is really necessary to 
diagonalise the full Hamiltonian of the N = L^ system to 
estimate the cost, A£, of a local move. Imagine a 'large' 
system, (with L = 20, say, for arguments sake) and some 
degree of 'disorder' seen by the electrons arising from 
the classical thermal fluctuations or quenched disorder. 
Quahtatively, if the effect of a change, r]i -^ ri[, does not 
'propagate' very far, as one would expect in a system 
with some disorder, then the energy cost of the move 
should be calculable from a Hamiltonian which involves 
only electronic degrees of freedom in the 'vicinity' of R^. 
We will discuss the analytic basis of such an argument 
separately, for the Monte Carlo it only requires that we 
modify step {ii) of the ED-MC strategy, discussed earlier. 
We compute Af as £c{'n'} ~ £c{'n}^ where £c{'ri} is the 
energy computed by constructing a Hamiltonian of Nc = 
L^ sites around R^, and diagonalising this Hamiltonian in 
the background configurations {ri}c and {ri'}c, where the 
curly brackets, {}c, refer to the configuration within the 



cluster. We will show results on 'equilibriation' within 
the TCA approach in the last section. 

After equilibriation the fermion properties are com- 
puted by diagonalising the full L^ Hamiltonian in the 
equilibrium background. Transport properties are calcu- 
lated based on the spectrum and states obtained from 
these diagonalisation, using the Kubo formula^^. 



III. RESULTS 



Apart from the electronic parameters, the two compu- 
tational parameters in the problem are N and Nc- In 
our notation TCA(iV : Nc) implies MC for a N site sys- 
tem based on a cluster of size Nc- ED-MC obviously 
corresponds to TCA{N : N). We assume a cube ge- 
ometry, with periodic boundary conditions (PBC) for 
both the system and the cluster. Ideally one should have 
TCA(A^ : N) available for large sizes and study conver- 
gence as Nc -^ N. Unfortunately ED-MC can be done, 
with great effort, only for sizes < 6"^, so TCA(iV : N) will 
be rarely available at large A'' and we have to analyse the 
approximation based on the following checks. 

(1). We study Hi and H2 using ED-MC on the largest 
possible lattice, A^ = 6"^. With somewhat reduced ther- 
mal averaging, and using scan in T (at fixed electron den- 
sity), or a scan in fj, (at fixed temperature), we establish 
the "large size" exact results. We then use TCA(A : Nc) 
with Nc — 3'^, 4"^ and 5'^ to check the convergence to 
TCA{N : N). 

(2). We study a "large" system, A^ = 8'^, and monitor 
the trend in TCA(A^ : Nc) with growing Nc, remaining 
in the regime Nc <^ N. We use Nc = 3^, 4"^, 5"^ to assess 
the convergence of the results to an asymptote with Nc 
still < N. 

(3). We compare the energy cost of actual microscopic 
MC updates between ED-MC and TCA. We evolve a 
system via TCA(8'^ : 4^) but simultaneously compute 
the energy cost of the updates via exact diagonalisation 
of the full 8'^ system. This is done by choosing a site 
randomly and computing both the exact 8'^ energy cost 
and the L^ energy cost whenever an update is attempted 
at that site. This yields information on how well TCA 
microscopically estimates the energy cost, rather than at 
the level of system averaged properties. 



A. Clean Holstein Model 



The Holstein model provides the minimal description 
of coupled electron and phonon degrees of freedom and, 
in the adiabatic limit, involves the following phases^^"^^, 
(i) a Fermi liquid (FL) metal, without any lattice distor- 
tions at T = 0, {a) a positionally disordered insulating 
polaron liquid (PL) at strong coupling, and (Hi) charge 
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FIG. 1. The result of fi variation in the clean Holstein 
model at coupling A = 2.5 and T = 0.02. The exact re- 
sult, with L = 6, is compared to results with Lc — 3 and 
Lc — 4. Panel (a), the variation in electron density, showing 
the discontinuous jump from the Fermi liquid to a charge or- 
dered state, (fe). the variation in charge order parameter with 
jj,, and (c) . the variance of the effective disorder (see text) seen 
by the electrons. We have not used Lc = 5 since the results 
on Lc = 4 are already very close to ED with 6^. 

ordered insulating (COI) phases close to n = 0.5. 
The physics of these phases has been discussed earlier 
within DMFT, and we also discuss it in detail in sepa- 
rate papers^*^'^^. Our intention here is to estimate the 
effectiveness of TCA in capturing the known features of 
the Holstein model as well as compare with exact MC 
calculations. 

Fig.l and Fig. 2 show the variation in carrier density, 
n(ii), the order parameter S'(7r,7r,7r) for commensurate 
charge ordering, and the variance of the 'effective disor- 
der', 77^ (defined further on), seen by the electrons, with 
varying fx at two temperatures. In a model which has the 
possibility of phase separation, and 'disallows' a certain 
density range, it is imperative to work with constant /i 
to map out the phase diagram. At the specified T and 
/i, TCA is used to obtain a family of equilibrium phonon 
configurations, which arc then used to solve the full elec- 
tron problem. 

Fig.l pertains to low temperature, T — 0.02, at inter- 
mediate coupling, where there are two phases, (i) a FL at 
low doping, n ^ 0.1, and (ii) a commensurate COI phase 
for 0.35 ^ n ^ 0.5, and a regime of phase separation for 
n between ~ 0.1 — 0.35. Before analysing the size depen- 
dence of the TCA results let us define the basic indica- 
tors. Fig.l. (a) shows n(/i), including the 'discontinuity' 
due to phase separation. Fig.l.(b) shows the COI order 
parameter S'(7r,7r, tt), computed from the structure fac- 
tor S{q) = {l/N^)E^A{n^){nJ))e''i■^^^-'^^\ where (n,) 



is the quantum average of n^ in a MC configuration and 
the outer angular brackets indicate average over config- 
urations. Fig.l.(c) shows the effective disorder based on 
the following prescription: for the Holstein model the 
electrons see a potential ^f ~ €i — Axf , where e^ is the 
extrinsic disorder and xf is the structural distortion in 
an equilibrium MC configuration [a is a MC configura- 
tion index). A crude measure of the 'disorder' seen by 
the electrons is provided by the variance of the ^i distri- 
bution, averaged spatially and over MC configurations. 
If we denote rji = ^i — ^, where ^ is the spatial average, 
then the effective disorder rj^ = (Vi)- ^ is the thermal 
and configuration averaged disorder that dictates the sin- 
gle particle scattering, and to some extent the transport 
properties. Our 77^ data in Fig.l.(c), and later figures, 
quantify this disorder. 

The data in Fig.l and Fig. 2 are on a 6^ system, using 
clusters of 3"^, 4^ and 6'^ itself to anneal the {xt}. Com- 
paring the cluster size dependence of the various physical 
quantities it is obvious that while updating using the 3^ 
cluster is unable to accurately capture the effects in the 
6'^ system, the results based on 4^ are qualitatively simi- 
lar to that of 6'^ . The size difference between the cluster 
and the full system is a factor of 6'^/4'^ r^ 3.4. The ra- 
tio of the computation time between ED-MC on 6'^ and 
TCA(6'^ : 4"^) is r^ 40, if the same extent of averaging is 
employed in both calculations. 

At higher temperature, T = 0.08, Fig. 2, the results 
based on 3^ continue to differ from the exact 6^ result 
but the 4^ result is virtually indistinguishable from the 





FIG. 2. Same as in Fig.l, except for T = 0.08. At 
this higher temperature the coexistence jump is almost 
smoothed out. The agreement between the exact result and 
TCA(6^ : 4^) is even better here. 
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FIG. 3. Holstein model in the metallic regime, A = 2.0 
and n = 0.3. Results with TCA(8^ : L^), with Lc = 3 - 5. 
Panel (a): resistivity p{T). Except at the lowest T, the results 
with Lc — 4: and Lc = 5 are virtually indistinguishable, while 
Lc — 3 matches the large cluster data only at high T. Panel 
(6): the effective disorder seen by the electrons, and panel 
(c) the density of states at the Fermi level. 

exact answer. The key to this lies in the large damping 
of the electrons arising from thermal fluctuations, and at 
this temperature the effect of the 4^ finite size gap is no 
longer relevant. 

Fig. 3 shows a different kind of result, where the sys- 
tem is studied via TCA at constant density, n = 0.3, at 
EP coupling A = 2.0, on a "large" system, L = 8, with 
cluster size varying from 3"^ — 5'^. We obviously cannot 
do an exact calculation on the 8'^ system, so the results 
in Fig. 3 are intended to check out (?) the convergence of 
the TCA data to the Nc — N asymptote with growing 
Nc, and (ii) study the temperature dependence of this 
convergence, since strong disorder, i.e, large 77^, at high 
temperature could make even the 3"^ based calculation 
viable. In Fig. 3 we directly compute the resistivity, us- 
ing a method described in an earlier paper^^, as well as 
rj^ and the density of states N{eF) at the Fermi level. 

Here again, the results based on 4*^ and 5"^ clusters are 
virtually indistinguishable except at the lowest tempera- 
ture. The result based on 3'^ is visibly different from that 
on 4'^ — 5'^ at the lower temperatures, but converges to 
a common value for T ^ 0.3, by which r/^ is quite large. 
The results on 77^ itself and N(ep) are quite similar for all 
Lc at all T, but the resistivity (which is a more stringent 
comparison) differentiates the changing character of the 
result with varying Lc- We would like to draw a general 
conclusion from these results, and back them up as we 
discuss the disordered Holstein model in the next section. 
If the single particle damping rate, F, arising out of 77^ 
is comparable to the finite size gap, 12t/L^, in the clus- 



ter then the specific finite size features of the cluster are 
smeared out and it mimics a 'large' system. So, annealing 
the variables on the large system via TCA is effective if 
r ^ W/Nc, where W is the bare bandwidth of the system. 
For a fixed T (and extrinsic disorder) this condition can 
be met by increasing Lc, while for a fixed Lc the accuracy 
increases as the net disorder (from thermal fiuctuations 
and extrinsic disorder) increases. 

To substantiate this claim, as well as check the ability 
of TCA to capture the 'fingerprints' associated with a 
specific disorder realisation, we next consider the disor- 
dered Holstein model. 



B. Disordered Holstein Model 



The key feature of the Holstein model is the possibilty 
of 'self-trapping', i.e, polaron formation, when the EP 
coupling exceeds a certain threshold. The critical cou- 
pling for single polaron formation -'^^ in 3d is Xc/t ^ 3.3, 
which implies that the polaron 'binding energy' E^ — 
\'^/{2K) ~ 5.44i. This would imply that at A = 2, where 
Ep = 2t we should be far from any polaronic instabil- 
ity. This is indeed true in the absence of disorder and 
Fig. 3, for example, shows that the response is metallic 
with dp/dT > 0. However, even weak disorder, A = 0.6, 
has dramatic effect in the FL phase, see Fig. 4. This figure 
presents results on the disordered Holstein model studied 
directly via ED-MC on 6^, as well as by TCA using 3^ 
and 4^ clusters on the 6^ system. 
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FIG. 4. Disordered Holstein: A = 2.0, A = 0.6, and 
n = 0.3, exact results using 6"^ and TCA. TCA is based on 
Lc — 3 and Lc = 4. Panel (a), resistivity p{T), panel (&). the 
effective disorder rj^ , panel (c). the density of states at the 
Fermi level. 
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FIG. 5. Thermally averaged density profile, top surface of a 
6'^ system: comparing exact results on 6^ with TCA. The top 
row is for Lc = 3, next for Lc — 4: and lowest for Lc = L — 6. 
Temperature, left to right, are 0.3, 0.11, 0.01. The disorder 
realisation {ei} is the same in all cases. 



In terms of the physical effect of disorder, the ED-MC 
indicates that the interplay of disorder and EP coupling 
can turn the system into a very bad metal (or even insu- 
lator) with a large resistivity at T = 0, and dp/dT < 
for T ^ 0, Fig. 4. (a). The 'effective disorder' seen by the 
electrons is large down to T = 0, Fig.4.(b), and there is a 
pseudogap in the DOS, as evident from N{ep), Fig.4.(c). 

The conversion of a FL (at A = 2.0, A = 0) into a 
'polaronic' phase by weak disorder happens because the 
density inhomogeneity created by weak disorder is dra- 
matically amplified by strong EP coupling^^'^^ leading 
to strong localisation. However, all electronic states are 
not strongly localised, as the spatial pattern. Fig. 5, and 
N{ep), Fig.4.(c), indicate. Fig. 5, discussed further on, 
shows the thermally averaged density n^. 

In contrast to our results on the clean Holstein model. 
Fig. 1-3, notice that all the sizes, 3^, 4^ and 6'^, yield sim- 
ilar results on all the indicators. Fig. 4. (a)- (c). While the 
4^^ based results almost coincide with the exact 6^ an- 
swer even the 3'^ based results capture all the qualitative 
features and even the numerical values reasonably accu- 
rately. The large effective disorder in the problem, arising 
from the strong lattice distortions, Xi, makes even the 3^ 
calculation acceptable. The rj^ in this model is ~ 3.0 
over the whole T range, comparable to the maximum rf 
in the clean problem, Fig.3.(b). 

While the transport and spectral properties seem to be 
adequately captured by TCA, does the method succeed 
in capturing the specific 'fingerprint' of a disorder reali- 
sation, {ei}? Fig. 5 shows the thermally averaged density 
pattern, n^, at three different temperatures (along the 
row) computed via TCA using Lc = 3 (first row) , Lc = 4 
(second row) and the exact, Lc = 6 case (bottom row). 
The TCA was run with the same realisation of disorder in 
all three cases. Remember that while the cluster based 
update is used for the phonon degrees of freedom, the 



final density field is calculated by diagonalising the full 
Hamiltonian in the background of the quenched disorder 
and the phonon configurations obtained via TCA. The 
cluster diagonalisation by itself cannot yield the density 
field. 

At first glance, the results of all three runs, compared 
along a column, match well. At intermediate and high 
temperature there is weak but still visible density con- 
trast and the results of all three schemes match very well. 
At the lowest temperature, the contrast is strongest and, 
although the correspondence along the third column is 
quite striking, there are some minor variations between 
the panels. This is partly because the MC based anneal- 
ing is less effective at low temperature due to the small 
acceptance rate of moves. Apart from this generic diffi- 
culty with MC calculations, we think the overall ability 
of TCA to capture the specific features of a disorder real- 
isation (and not just system averaged properties) is quite 
impressive. 

Now consider using TCA on large sizes, L = 8, as 
in Fig. 3, for the disordered problem. Fig. 6 shows TCA 
based results using Lc = 3, 4, 5 to solve the L = 8 prob- 
lem, and Fig. 7 shows the associated density profile. As 
stated before it is impossible to do ED-MC on 8"^, so 
the data here is meant to indicate the convergence of the 
TCA based results to the Nc ^ N asymptote, although 
the exact result at A^c = A^ is not available. The effective 
disorder, Fig.6.(b), and N{eF), Fig.6.(c), are very simi- 
lar for all Lc, and even the resistivity. Fig. 6. (a), matches 
quite well after disorder average over ^4 — 5 copies. 
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FIG. 6. Disordered Holstein model, A = 2.0 and n — 0.3, 
A = 0.6. Results witli TCA(8^ : L'i), witli L^ = 3 - 5. Panel 
(a): resistivity p{T), (6): the effective disorder seen by the 
electrons, and (c) the density of states at the Fermi level. 
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FIG. 7. Thermally averaged density profile: comparing 
TCA based results on a 8"^ system using Lc = 3, 4, 5. The top 
row is for Lc — 3, next for Lc = 4 and lowest for Lc — L = 5. 
Temperature, left to right, are 0.4, 0.14, 0.04. The disorder 
realisation {e;} is the same in all cases. 

We would argue that for the conditions studied, TCA 
with Lc ~ 3 — 5 is quite adequate in annealing the phonon 
variables on i = 8. In fact when Lc is large enough, 
so that F J^ 12t/L'^, the outer limit, i.e, system size L, 
is actually irrelevant. Updates based on small clusters 
can successfully generate the appropriate configurations 
on large lattices. The only reason large L is needed at 
all, for the system as a whole, is to compute transport 
properties, or check for long range order. 

The thermally averaged density profile corresponding 
to the MC in Fig. 6 are shown in Fig. 7. As in Fig. 5, 
the patterns along each row correspond to decreasing T, 
reducing from T = 0.4 to 0.14 to 0.04. The first row 
corresponds to TCA with Lc = 3, the second with Lc = 4 
and the third to ic = 5. All three systems have the same 
realisation of quenched disorder {e^}. Again, as in Fig. 5, 
the different TCA results are in excellent agreement at 
higher temperature (first two columns) while there are 
some differences at the lowest temperature due to the 
difficulty in annealing. The actual percent difference in 
the density field, between the three systems, averaged 
over the whole 8"^ system is ~ 10 — 15%. 



C. Double Exchange Model 

The clean double exchange model has been widely 
studied, using a variety of analytical approximations and 
numerical techniques. We do not enter into a detailed 
recapitulation of these results since we have already re- 
viewed them in detail^^ in an earlier paper. Here we focus 
primarily on MC based results, since these are unbiased, 
although necessarily finite size. The ground state of the 
clean DE model is a saturated ferromagnet and in terms 



of transport both the ferromagnetic and paramagnetic 
phase are metalhc. 

In the hmit Jn/t ^ oo the DE model maps on to a 
spinless fermion problem, with 'hopping disorder' arising 
from the background spin configuration. Since the strong 
local coupling Ju couples the core spin orientation to 
fermion spin projection, only the 'locally aligned' fermion 
state is viable at each site (the other is at an energy Jh 
above). The electron hopping between two sites depends 
on the electronic eigenfunctions at each site, and so on the 
spin orientation, 
out to be: 



The 'projected' Hamiltonian turns 



H ^-t^{ 9^J-^l-^J +h.c)- ii^n, 

= -i ^ /,, ( e^*'^ 7h. +h.c)-^iY,n, (2) 



(»j) 



The 7's are spinless fermion operators. The hopping 



amplitude gij = fi^ 



between locally aligned states, 



can be written in terms of the polar angle (6i) and 
azimuthal angle (0i) of the spin S^ as, cos-j-cos-i- 
+sin^sin-j-e~^ (0*-"^^). It is easily checked that the 
'magnitude' of the overlap, /^ — ^{1 + Si.Sj)/2, while 
the phase is specified by tan^ij ~ Inn{gij) / Re{gij) . 

Fig. 8 shows MC results on the DE model, using ic = 4 
and varying the system size from L — A to L ^ 10. This 
is unlike our earlier results on the Holstein model where 
we kept the system size fixed and varied Lc, looking for 
convergence. For the DE model, there is already MC 
data available^^ on large sizes. Fig. 8 shows the evolution 
of the magnetisation profile m{T) as L increases. 




FIG. 8. Magnetisation in the clean DE model with Lc = 4, 
and L = 4 — 10. The electron density is n = 0.32. The 
magnetisation profile obtained by extrapolation to 1/ ^ oo 
is shown as a dotted line. The arrow indicates Tc using the 
the moment expansion based MC and finite size scaling by 
Furukawa et al.. The inset compares the TCA result to an 
earlier 'effective Hamiltonian' approximation made by us. 
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FIG. 9. Disordered DE model: magnetisation and resistiv- 
ity, with Lc — 3, 4, 5, and L = 8. The electron density is 
n = 0.30. 

The large L extrapolation of this trend, at Lc = 4, 
suggests Tc w 0.14. The result of moment expansion 
based MC and finite size scaling^^ indicate Tc « 0.12 
at this density. Although there is a difference between 
these results, the ability of TCA to approach the exact 
answer, with substantially less computational effort, is 
evident. The inset in Fig. 8 compares an earlier approxi- 
mation used by us^'"' (using a classical effective Hamilto- 
nian for the spins) with the TCA using L = 10, Lc ~ 4. 
The improvement in TCA, particularly in capturing the 
Tc scale is obvious. 

Fig. 9 shows results on the DE model in the presence 
of weak disorder in the electron system. This model has 
an additional term ^^ e^ni where the e^ is binary disor- 
der with strength ±A as in the Holstein problem. The 
data on magnetisation and resistivity in Fig. 9 is shown 
for N = 8"^, with Lc ~ 3,4,5. As in the Holstein prob- 
lem, the profiles with Lc — A and Lc ~ 5 are barely 
distinguishable. 



D. Computational indicators 

L Equilibriation 

Since the TCA based MC does not compute the total 
energy of the system in the process of updating, unlike 
ED-MC, the process of equilibriation and stability of the 
energy is not obvious. To explicitly confirm the nature of 
equilibrium fluctuations (and the absence of drift in the 
mean value) as well as visualise the process of equilibria- 
tion in response to a temperature step, we ran TCA with 
a simultaneous calculation of the total energy at the end 
of each MC sweep. The results, for the Holstein model, 
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FIG. 10. Total energy variations during the TCA-MC on 
the Holstein model. The electronic parameters and temper- 
atures are marked on the Figure. This trace is quite similar 
to that obtained in a typical ED-MC run, and the energy 
fluctuations reduce as expected with decreasing temperature. 
The inset shows the thermally averaged energy at each T as 
a function of temperature. 

are shown in Fig. 10. Notice that this requires diagonal- 
isation of the full Hamiltonian matrix Nmc times, rather 
than NmcN times as required by ED-MC. Apart from 
assurance about proper equilibriation this diagnostic al- 
lows us to track potential hysteresis effects, i.e, difference 
between heating and cooling, if a first order transition is 
involved. 



2. Comparing A£ 

It may be useful to check out the claim that microscop- 
ically the TCA based calculation of A£ quickly converges 
with increasing Nc (even if it is still <^ N) and with in- 
creasing extrinsic disorder. To that effect Fig. 11 analyses 
the results on the Holstein model. We update the sys- 
tem using TCA with Lc = 3 — 5onaL = 8 system. 
For each Lc, and specified T, we choose a reference site 
randomly, and whenever a phonon update is attempted 
at that site in the course of the MC sweep, we not only 
compute the cluster based energy difference A£c but also 
the exact energy cost of such a move A£ed based on the 
full system. 

The system evolves accoring to TCA, but we keep track 
of these energy differences (at that site) and construct the 
following error measure: 



S{L,,L) 



1 , Af , - ASed , 



N 



MC 



AS 



ED 



(3) 



where the averaging is over MC steps, and the 'error' 5 
implicitly depends on temperature, as well as all other 
electronic parameters (in particular, disorder). 

Fig. 11 shows results, again on the Holstein model, with 
L = 8 and Lc — 4, panel (a). The state at T = is a clean 
Fermi liquid for T = 0, the electronic states are simple 
tight-binding states and, as expected, AS computed on 
8'^ and 4"^ have a fair difference, about 25%. This "error" 
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FIG. 11. Error in energy estimate of MC moves. The indi- 
cator S defined in the text, compares energy costs on the full 
L^ system with that based on L^ clusters. Panel (a) shows 
the T dependence of 5 with increasing disorder, and the inset 
shows the Lc dependence at A = 0. Panel (6) same as in 
panel (a) with different parameters, marked on the Figure. 
Inset to panel (b) variation of 5 with the effective disorder, 



"^ combining the results in main panel (a). 



falls quickly with increasing temperature, as thermal fluc- 
tuations in the Xi increase (staying at A = 0), and also re- 
duces systematically on introduction of even weak disor- 
der A = 0.3 — 0.6. The inset to Fig. 11. (a) shows temper- 
ature dependence of the error with varying Lc at A = 0. 
To construct an approximate measure, the typical error, 
averaged over the temperature window is ~ 0.5 at Lc = 3, 
-- 0.1 at Lc = 4, and ^ 0.05 at Lc = 5. Fig.ll.(6) shows 
S a.t n — 0.5, A — 2.0 (which is a charge ordered state at 
low temperature) and n = 0.3, A = 3.0 (which is a pola- 
ronic insulator) again at L = 8 with Lc = 4. The error 
in the CO problem is non monotonic because, unlike the 
Fermi liquid, the system goes into a ordered chessboard 
phase at T = and this localisation reduces the error. 
In the polaronic insulator phase the error is below the 
5% threshold at all temperatures due to the strongly lo- 
calised nature of electronic wavefunctions. The inset to 
Fig. 11. (6) puts together the error variation with respect 
to the effective disorder, 77^, for the data in Fig. 11. (a) 
for A = — 0.6 with varying temperature. As we have 
argued earlier there is a roughly 'universal' behaviour of 
the error in terms of the effective disorder, irrespective 
of its origin,' and at Lc = 4 an error < 0.1 is obtained 
whenever 77^ > 1.0 

Note that S above is a measure of error in the energy 
estimate, we have checked that the error in e~^^^ itself, 
which controls the MC moves, between the 4*^ and 8^ 
cluster are much smaller than the error in AS. This 
explains why the thermodynamic and transport results 
are accurate even with Lc = 4 in the clean system. 



IV. CONCLUSION 



Since we have put forward the TCA based approach to 
adiabatic problems as a 'general' many body technique, 
not restricted to the specific examples studied here, let 
us compare it with DMFT^ which finds wide use as a 
method for handling correlated electrons, (i) DMFT 
maps on a correlated lattice problem to an impurity 
model, reducing it to an effective single site problem and 
making it more tractable. The adiabatic approach, and 
the TCA approximation, ignores the quantum dynamics 
of the background fields, leading to a static but 'annealed' 
disorder problem, (m) The two methods have comple- 
mentary reach. DMFT is good at handling local quan- 
tum correlations but misses out on spatial fluctuations 
and disorder physics. The real space adiabatic approach, 
supplemented by TCA like approximations, handles dis- 
order and spatial correlations accurately but cannot han- 
dle interacting quantum modes, (in) Both have had sig- 
nificant success in materials physics. DMFT extended to 
include non local correlations will be truly global method, 
the adiabatic approach extended to include quantum dy- 
namics, and implemented via TCA like approximations, 
or refinements, would 'reach' real life problems from an- 
other direction. 

In conclusion, we have put forward and benchmarked 
a new Monte Carlo technique for handling fermions 
strongly coupled to classical degrees of freedom. The 
method relies on the approximate 'locality' of the energy 
cost of a MC move in a system with moderate disor- 
der, allowing accurate estimate of energy differences to 
be made using on a cluster Hamiltonian instead of the full 
system. This allows annealing of the classical variables 
on large lattices, breaking the A''^ barrier that plagues 
exact diagonalisation based Monte Carlo. The approach 
makes no assumptions regarding the the starting Hamil- 
tonian, except the quadratic nature of the quantum de- 
grees of freedom. We have provided basic benchmarks 
on the Holstein and Double Exchange model, and will 
provide several new results in forthcoming papers. 

We acknowledge use of the Beowulf cluster at HRI. 
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